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Abstract 

In the last few years the realization has emerged that the universal baryons are 
almost equally distributed by mass in three components: (1) galactic concentra- 
tions, (2) a warm-hot intergalactic medium (WHIM) and (3) a diffuse intergalactic 
medium. These three components are predicted by hydro dynamical simulations and 
are probed by QSO absorption lines. To observe the WHIM in neutral hydrogen, 
observations are needed which are deeper than log(Nf//)=18. The WHIM should 
appear as a Cosmic Web, underlying the galaxies with higher column densities. We 
have used the WSRT, to simulate a filled aperture by observing at very high hour 
angles, to reach very high column density sensitivity. To achieve even higher image 
fidelity, an accurate model of the WSRT primary beam was developed. This will 
be used in the joint deconvolution of the observations. To get a good overview of 
the distribution and kinematics of the Cosmic Web, a deep survey of 1500 square 
degrees of sky was undertaken, containing the galaxy filament extending between 
the Local Group and the Virgo Cluster. The auto-correlation data has been reduced 
and has an RMS of ANhi = 4.2 x 10 16 cm -2 over 20 kms -1 . Several sources have 
been tentatively detected, which were previously unknown, as well as an indication 
for diffuse intergalactic filaments. 
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1 Introduction 



The number of detected baryons in the Low Redshift Universe is significantly 
below expectations. According to cosmological measurements the baryon frac- 
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tion is about 4% at z ~ 2 (Bennett et all (120031 ) : ISpereel et al . (2003) ). This is 



consis te nt with actua l numbers of baryons detected at z > 2 ([Weinberg et al. 
(119971 ); iRauchl (119981 )). In the curren t epoch howeve r , at z ~ about half 



of this matter has not bee n observed (Fukugi ta et al. (1998) 



(Il999h : iTripp et al.l (l2000h : ISavage et all (120021 ): IPenton et all fo004 )) 



Cen fc Ostriker 



Recent hydrodynami cal simulations g ive a possible solu t ion f o r the Miss- 



i ng B aryon problem (ICen fc Ostrikerl (119991 ); IDave et al.l (120011 ); iFang et al. 



(120021 )). These simulations show that in the current epoch, cosmic ba ryons 
are almost equally distributed amongst three phases (IDave et al.l Il999l ): (1) 
diffuse phase, (2) a shocked phase and (3) a condensed phase. The diffuse 
phase is associated with warm, low-density photo-ionized gas. The shocked 
phase consists of gas that has been heated by shocks during structure for- 
mation, with a moderate overdensity. Because the temperature range is very 
broad from 10 5 to 10 7 K, it is also called the Warm Hot Intergalactic Medium 
(WHIM). The condensed phase is associated with cool galactic concentrations. 
These three components are each coupled to a decreasing range of baryonic 
over-density: log(pH / Ph) < 1, 1-3.5, and > 3.5 and are probed by QSO ab- 
sorption lines wi th specific ranges of ne utral column density: log (iV#j) < 14, 
14-18 and > 18 dBraun fc Thilkerll2005h . 



2 WHIM 



The Warm Hot Intergalactic Medium is thought to be formed during struc- 
ture formation. Low density gas is heated by shocks during its infall onto the 
filaments that define the Large Scale Structure of the Universe. Most of these 
baryons are still concentrated in unvirialized filamentary structures of highly 
ionized gas. 



The W HIM has been obser vationally detected in O VI absorption (ITripp et al 
2004 ). but also via Ne VIII (ISavage et al.ll2005l ) and X-ray absorption (INicastro et al 
20051 ) . Of course absorption studies alone, do not give us complete information 
on the spatial distribution of the WHIM. Emission from the WHIM would give 
entirely new information about the distribution and kinematics. 



Direct detection of th e WHIM is very difficult in the EUV and X-ray bands 
(ICen fc Ostrikerl Il999l ). The gas is ionized to such a degree, that it becomes 
"invisible" in infr ared, optical or UV light, but should be visible in the FUV 
and X-ray bands (INicastro et al.ll2005l ). Given the very low density, extremely 
high sensitivity and a large field of view is needed to ima ge the filaments- 



Capab l e detectors are not yet available for the X-ray or FUV (lYoshikawa et al. 
(l2003h : INicastro et all (|2005h ). 
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Fig. 1. Observing mode of the WSRT dishes, a filled aperture of 300 m can be 
simulated when observing at large hour angles. 



Large -s cale gas filaments h a ve been detected in X-ray e mission (iWang et al 
f ll997h : IScharf et all fl2000h : iTittlev fc Henriksenl tOOlh ). With the current 
facilities, these X-ray emission studies reveal gas which is hotter and denser 
than the WHIM. This gas is not exp ected to contain a substantial part of 
WHIM baryons at the current epoch (jDave et all 1200 ll ). The WHIM should 
be most prominent in the soft X-ray band. The interstellar medium of the 
Galaxy complicates its detection since in the soft X-ray band the Galaxy is a 
strong source of emission at an effective temperature of ~ 10 6 K. The neutral 
hydrogen in the Galaxy on the other hand, prevents observations of the EUV. 



We have adopted a different approach by observing the WHIM in neutral 
hydrogen emission. 



At the current epoch we can confidently predict that in going down from HI 
column densities of 10 19 cm -2 (which define the current edges of well studied 
nearby galaxie s in HI emission) to 10 1 7 cm ~ 2 the surface area will increase by 



a factor of 30 (jBraun Sz Thilkerll200l . The critical observational challenge is 
crossing the "HI desert", the range of log(N^/) from about 19.5 down to 18 
over which photo-ionization by the intergalactic radiation field produces an 
exponential decline in the neutral fraction from essentially unity down to a few 
percent (eg. iDove fc Shuhl fll994f )h Nature is kinder again to the HI observer 
below log(N^/) = 18, where the neutral fraction decreases only very slowly 
with log(N^/). The neutral fraction of hydrogen is thought to decrease with 
decreasing column de nsity from about 10 0% for log(N HI ) = 19.5 to about 
1% at log{N HI ) = 17 dDove fc ShullHl994T ). The baryonic mass traced by this 
gas (with a 1% or less neutral fraction ) is expected to be comparable to that 
within the galaxies, as noted above. 
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degrees 



Fig. 2. The Point Spread Function (PSF) of the primary beam of the WSRT on 
logarithmic scale, averaged over a full band. 



3 Observations 



To obtain the highest possible brightness sensitivity, the WSRT has been 
used to simulate a large filled aperture. 12 of the 14 WSRT 25 m telescopes 
are positioned at regular intervals of 144 m. When observing at very low 
declinations and extreme hour angles, a filled aperture can be formed (as can 
be seen in Fig. 1), which is 300 x 25 m in projection. In this peculiar observing 
mode the excellent spectral baseline and PSF properties of the interferometer 
are still obtained while achieving excellent brightness sensitivity. A deep fully- 
sampled survey of the galaxy filament joining the Local Group to the Virgo 
Cluster has been undertaken, extending from 8 to 17 hours in RA and from -1 
to +10 in Dec. and covering 40 MHz of bandwidth with 8 km s _1 resolution. 
Mosaic-mode drift-scan observations are undertaken twice for each of 22,000 
positions, once at positive hour angle and once at negative hour angle. The 
expected survey RMS sensitivity is AN HI ~ 2 x 10 17 cm~ 2 over 20 km s _1 in 
the synthesis data. 



3.1 Primary Beam 



To obtain high dynamic range and image fidelity in a mosaiced field, the shape 
of the primary beam must be known to high precision. The approximation for 
the beam shape used in the past has strong limitations, since it assumes the 
beam to be axisymmetric and is truncated at a level of a few percent, which 
is above the level of the first sidelobes. We modeled the WSRT primary beam 
in great detail. Three series of holographic measurements were undertaken to 
fully sample angular scales between a few arcminutes and several degrees. An 
empirical model-cube was constructed, containing the primary beam pattern 
between 1322 and 1457 MHz, with a 1 MHz frequency resolution. In Fig. 2 
the model is plotted on a logarithmic scale, averaged between 1424 and 1440 
MHz. In Fig. 3 (left panel) the integrated main beam and side lobes are plotted 
on a normalized scale against frequency. There is a very significant periodic 
modulation of the beam size of ~ 4% in addition to the expected linear scaling 
with frequency. This modulation has a period of ~17 MHz. When doing the 
joint-deconvolution of the mosaic observations, an interpolated PSF will be 
calculated, corresponding to the frequency of the data. 

Compared to the previously employed cos 6 function, the new model has a sys- 
tematically broader main beam (at all phases of the oscillations), a "diamond" - 
shaped departure from circular symmetry as well as inclusion of the first few 
near-in sidelobes (and their strong variation with frequency). When correcting 
for the primary beam attenuation, this can have a significant influence on the 
calculated source fluxes, especially far from the beam center. To investigate 
the difference between the new model and the old beam description we have 
compared the implied flux densities of off-axis continuum sources with those 
obtained from the NVSS survey. 



The NVSS (NRAO VLA Sky Survey: ICondon et al.l (Il998h ) is a large area 



radio survey at 1.4 GHz. The images are on-line an d can be download ed. We 



have employed the WSRT SINGS continuum survey (IBraun et al.ll2006l ) which 
contains single pointing observations of more than 30 fields at the same fre- 
quency which means that they can be easily compared. We have produced an 
inital sample of 36 compact sources for which fluxes have been determined in 
both the NVSS images as well as the WSRT single pointing observations. The 
WSRT observations have been corrected for primary beam attenuation with 
the old cos 6 approximation and with the new model. The sources are located 
at different distances from the pointing center, to determine the reconstructed 
flux accuracy at different radii. 

The results are plotted in Fig. 3 (right panel), where the obtained fluxes are 
plotted against the radius from the center of the beam. All fluxes are nor- 
malized to the fluxes of the NVSS observations which are used as a reference. 
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Fig. 3. Left: The ~17 MHz period for the integrated power of the main lobe and 
the side lobes. Right: The normalized NVSS fluxes with the old and new beam 
descriptions, especially at large radii from the beam center, the new model gives a 
significantly better correction for the Primary Beam Attenuation 

At small radii, the old approximation (stars) as well as the new model (open 
circles) are very comparable with the NVSS-fluxes (straight line). When going 
to larger radii, there are more differences. The main beam of the cos 6 approx- 
imation is clearly too small, which means that the correction for attenuation 
will be larger and flux discrepancies will be larger when going out from the 
center. Indeed, at large radii the old approximation gives fluxes which are 
significantly too high, while the new model agrees with the reference values. 
Much larger samples of comparison sources will be required to document the 
accuracy of the new correction in detail. 



4 Results 



All observations for our WSRT Virgo Filament Survey (WVFS) have now been 
completed. We have begun analysis of the database by first reducing the auto- 
correlations which were obtained simultaneously with the cross-correlation 
mosaic. The first data cubes and moment images are now available for this 
total power data. Figure 4 shows the integrated column density map of the 
observed region. In addition to the many hundreds of known galaxies which 
were detected in this region there are a total of nine discrete sources found, 
which have not been detected before. All sources were selected to be above a 
level of 7a. The radial velocities are corrected for Virgo infall and the distances 
are calculated, using H = 71 km s _1 Mpc -1 . Global properties for the nine 
objects are given in table 1. 

Between the galaxies, diffuse filaments are apparent with column densities 
going down to ~ 8 x 10 17 cm -2 . While this is exactly the type of signature 
we were hoping to detect, further analysis will be required to rule out an 
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Global properties of the detected sources. 
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Fig. 4. Column density map with contour levels at 8 x 10, 3 x 10 18 and 3 x 10 
cm' 2 . The new objects are encircled and assigned a letter. 

instrumental origin for these features. 



5 Discussion 

The Warm Hot Intergalactic medium is expected to contribute about 30-50 % 
to the total baryon budget. Observing the WHIM is very difficult and ex- 
tremely sensitive methods are needed. We are observing the WHIM in neutral 
hydrogen emission since the expected column densities, between log(N#/) 14- 
18, are just accessible to current instrumentation at the high end. The first 
results look promising, since we reach a sensitivity of ANhi = 4.2 x 10 16 cm -2 
over 20 km s -1 in the auto-correlation data of our survey. 
Nine discrete sources have been tentatively detected which were previously 
unknown, as well as an indication for diffuse filaments that may be the HI 
counterparts of the " Cosmic Web" . 

Once the distribution of the WHIM is known, absorption lines toward strong 
background sources will be sought. These may provide both the metalicity and 
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the neutral fraction. With the neutral fraction in hand, the total mass of the 
WHIM can be estimated, allowing us to test the conjecture that the "Missing 
Baryons" have finally been found. 
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